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1. Introduction 



1.1 Background 

Because of the environment in which they operate, the 
seakeeping characteristics of a ship are of paramount impor- 
tance when assessing its overall performance. In the past 
this aspect of a ship's performance had to be judged by the 
results of model tests conducted at a point in time well into 
the preliminary design phase. While these tests provide 
results of good quality, they were not obtained until the 
pending design was quite firmly established. In fact, the 
results obtained by model tests had the characteristic of 
being just that, results, rather than an important input into 
the design cycle. The obvious desire then would be to have 
a tool capable of providing accurate predictions of seakeeping 
performance based on the data available in the conceptual 
design phase. These predictions could then be used to in- 
fluence the selection of hull form coefficients, etc. prior 
to the time when the hull form is actually being generated. 

In 1975 Professors T. Loukakis and C. Chryssostomidis 
published the "Seakeeping Standard Series for Cruiser-Stern 
Ships". ^ This paper corrolates the seakeeping behavior, 
as predicted by computer model, of the Extended Series 60 
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hull forms and sets forth a method by which the performance 
of this type of ship may be predicted based on five para- 
meters: Froude number, F; ratio of significant wave height 

to ship length, S; beam/draft ratio, B/T; length/beam ratio, 
L/B? and block coefficient, CB. With this procedure a 
designer can predict the relative merits of various candidate 
designs at a very early stage. This represents a significant 
capability. 

As a result of the work represented in reference [1] 
there is considerable interest in generating a similar sea- 
keeping series for contemporary cruiser /destroyer type hull 
forms. In order to do this in the fashion of reference [1], 
a representative sample of the ship type must be analyzed by 
computer model and then the results corrolated. It was 
this need for sample hull forms that provided the motivation 
for this thesis. 

1.2 Thesis Content 

There are two aspects of hull form generation 
addressed in this thesis: first, hull form modification and 

second, mathematical lines representation and fairing. The 
technique of hull form modification developed in chapter two 
is based on the work of H. Lackenby reported in reference [2] 
The essence of this method is that the sectional area curve 
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of an existing ship is redrawn in a systematic fashion to 
produce a curve with the desired values of prismatic co- 
efficient, C p , and longitudinal center of buoyancy, LCB. The 
sections are then shifted longitudinally to produce a modified 
form with these characteristics. In applying this technique 
to destroyer type ships there were several anomalies encoun- 
tered which required that the method of Lackenby be further 
modified. These modifications, with the pertinent background 
are contained in chapter two. 

The other aspect of lines generation which is addressed 
in chapters three and four is mathematical lines representation 
and fairing. Although the fairness of a hull form is not 
critical to the seakeeping analysis it is an unavoidable 
subject when considering computer aided ship design. In 
these chapters the use of parametric cubic splines and least 
squares curve fitting are addressed. While the parametric 
splines are shown to provide the capability of representing 
virtually any type of line, the least squares fairing technique 
is limited to use with curves representable by single valued 
functions. The algorithms are, however, capable of fairing 
lines with infinite slopes at the end points. 

It is anticipated that the tools developed in this 
thesis could be readily fused into a single computer program 
with the capability of modifying an existing ship form to 
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obtain a faired design with the desired coefficients of form. 
When this is developed it will be possible to generate 
rapidly any number of designs for subsequent performance 
analysis. The implications of this are discussed in chapter 
six. 
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2. Method of Hull Form Modification 



2.1 Background 

During the design of all but the most trivial engin- 
eering systems, it is incumbent upon the engineer that he 
or she formulate a model of that system. Additionally, 
the designer must continually refine the model with each 
successive iterative cycle so that the results are of 
sufficient detail to be meaningful [ 3 ] . One such model 
used during ship design is a geometric description of the 
ship's hull form. The most traditional manner of providing 
this information is by way of the lines drawing. 

The ship's lines drawing, more frequently referred 
to as the ship's lines, is a set of three orthogonal views 
of the ship's hull depicting the lines of intersection of 
various planes with the hull form. When viewed in conjunc- 
tion with one another, they provide the capability to 
spatially locate any point on the moulded surface of the 
ship. Figure 2.1, taken from reference C 4 ] #> is an example 
of a lines drawing for a "Mariner"-class, steel hull cargo 
vessel . 

While the lines drawing, prepared manually by the naval 
architect and draftsman, has been the older and more tradi- 
tional means of depicting the geometric properties of a ship. 
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the advent of the high speed computer has provided great 
impetus for defining the ship's form in a mathematical 
format [5]. It is interesting however , that a very 
successful attempt was made at representing ship's lines 
mathematically by Admiral David Taylor in the early 1900 's 
[ 6 ] . This will be expanded upon a bit later. 

When first confronted with the job of creating the 
lines of a new ship, the naval architect seeks a means of 
quantifying the expected form of the vessel so that he may 
strive to create an "optimum'' design. These optimizing 
criteria generally take the form of requirements and restric- 
tions placed on the various coefficients of form, i.e., 

C , C , LCB, LCF, etc. However, there might also be 
p w 

requirements placed on certain specific regions of the 
ship. An example of this could be the shape of the midships 
section for a cargo vessel or the stern configuration 
dictated by propeller and rudder selections. Nonetheless, 
when the naval architect completes his candidate design, 
the important product will be a faired set of ship's lines 
meeting all the optimizing criteria previously established. 

The above procedure is clearly long and involved. For 
this reason much effort has been expended to develop hull 
form modification techniques. The objective of these pro- 
cedures is to utilize an existing, successful hull design, 
or parent form, as a basis and then to alter this form in a 
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systematic fashion. This modified form should have the 
desired characteristics, and, hopefully, require little 
additional refairing. It is this fairing procedure, des- 
cribed in chapter 4, which requires a large proportion of 
the designer's effort. The remainder of this chapter 
addresses the modification techniques themselves. 

One of the oldest and most widely used methods of hull 
form modification is illustrated in Figure 2.2. In essence, 
the sectional area curve of an existing ship is altered by 
some arbitrary or systematic method to produce a curve which 
satisfies some criteria of the designer, usually prismatic 
or block coefficient and longitudinal center of buoyancy. 

The offsets for the new design are then obtained by taking 
the section in the parent whose ordinate in the sectional 
area curve matches the ordinate of the derived curve. This 
is represented by the movement of section A at position X a 
in the parent to section A' at position X a , in the derived 
hull form. Hence, it is merely a longitudinal repositioning 
of the existing sections. This method works reasonably well 
as long as the shifts are of "moderate" amount and the 
designer is prepared to accept the resulting profile and 
water lines without alteration. 
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The above method lends itself quite well to design 
without the aid of computers or other automatic computational 
devices. However, in recent years there has been much work 
done in the area of hull form modification with the use of 
high speed digital computers. In virtually all cases where 
computers are used, an effort is made to represent the ship's 
contours or surface regions [ 5 ] in a mathematical format. 

It is for this reason that the "Taylor Standard Series" is of 
interest. It was Admiral David Taylor who, in the early 
1900's, generated one of the first successful hull series 
based on representing the sectional area curve and design 
water lines by fifth degree polynomials [ 6 ] . 

Another procedure of hull form modification utilizes 
specific transformation functions to alter various regions 
or characteristics of the hull [ 7 ] . This form of 
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modification provides the user with much greater control over 
the specific ship form than the method of shifting sections 
longitudinally as previously described. An interesting 
description of this type of procedure may be found in 
reference [ 7 ] . 

It is, however, the method of longitudinally shifting 
sections which was selected for development in this thesis. 

The reasons for its selection are twofold. First, preliminary 
work with destroyer- type ships conducted at M.I.T. during the 
summer of 1977, indicated that the results of the modification 
were quite realistic and not plagued by gross unfairness. 
Secondly, the method was tractable and readily adopted to the 
peculiarities of destroyers. Those peculiarities being 
principly the fact that this type of ship has no parallel 
middle body and also that the section of maximum area, in 
most cases lies at a location other than midship. 

The specific method of modification used is that of 
Lackenby [ 2 ] as subsequently modified by Moor [ 8 ] , and then 
again by this author. Briefly, the developments presented in 
reference [ 2 ] are highly general, permitting the designer to 
vary the value of prismatic coefficient, C^, and the longi- 
tudinal center of buoyancy, LCB, of a very wide variety of 
ships. Included was the capability of altering, or retaining 
unchanged, the parallel middle bodies of ships so configured. 
However, one serious drawback was that the method left no 
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control over the design waterline, and while this line might 
turn out fair, the longitudinal center of floatation, LCF, 
merely ended up where it did. It was to this problem that 
Moor [ 8 ] was concerned. By addressing himself to both the 
sectional area curve and the design waterline in the manner 
of Lackenby, and then coupling the two procedures, he was 
able to obtain a derived form having the desired values of 
C , C w , LCB and LCF. 

At this point only one minor problem existed with the 

method as it stood. For ships with keelrise fore or aft it - 

was possible to obtain unwanted oscillations in the ship's 

centerline profile. To eliminate these oscillations, this 

author extended the logic of Moor to include the ship's 

profile. In so doing, the designer may be assured of a 

derived form having not only the four desired characteristics 

and coefficients previously mentioned, but also the desired 

profile. The only restricting requirement, other than the 

fact that the changes in C and C be "moderate", is that 
p w 

for the method to be mathematically rigorous the section 
of maximum beam, sectional area and local draft must be 
coincident. If this isn't the case a small (*1%) unpre- 
dictable error, based on the parent hull design and the 
desired changes is introduced. 
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All of these relationships are developed in full in 
the following section. The only other alteration to the 
methods of Lackenby and Moor was that the procedure had to 
be capable of accommodating destroyer-type ships whose 
maximum sections fill other than at midship. This change 
is also included in the derivations that follow. 
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2.2 Development 



2.2.1 The "One Minus Prismatic" Variation 

As a means of introducing the method of longitudinal 
repositioning of sections, the traditional "one minus 
prismatic" is first developed. This procedure enables the 
designer to modify the fineness of a parent form by expanding 
(creating in ships without) , or reducing the region of 
parallel middle body. It is convenient for this, and the 
following derivations to refer to Figure 2.3 and the 
following definitions. It should also be noted that the 
sectional area curve is normalized with respect to both the 
value of maximum area and length of the half body. 




Figure 2.3- S.A. C 
(with parallel middlebody) 



Figure 2.3- S. A. Curve 
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For the parent design: 

<J> = the prismatic coefficient of the half body, 
x = the fractional distance from midships of 
the centroids of the half body, 
p = the fractional parallel middle of the half body, 

x = the fractional distance of any transverse 

section from midships. 

y = normalized area of any transverse section at 
longitudinal position x. 

For the derived form: 

5<j> = the required change in prismatic coefficient 
of the half body. 

5p = the resulting change in parallel middle body. 

fix = the necessary longitudinal shift of the 

section at x required to generate the required 
change in prismatic coefficient, 
h = the fractional distance from midships of the 
centroid of the added "sliver" of area 
represented by 5$. 

In Figure 2.3 it should be recognized that AB'C is the 
curve of the derived form and curve ABC is that of the 
parent. In accordance with the method of the "one minus 
prismatic" the new location of the transverse sections is 
defined by the following equations. 
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1- (X+ 6x) _ 1- (0+60 ) 
1-x “ 1— 4> 



6x _ 60 
I=x “ I=* 



ax ■ ( i_x > 



( 2 . 1 ) 



The area BCD is seen to be 1 - 0 and the above modification 



The new area B*CD is therefore 1 - (0+60), demonstrating 
that the method generates the desired prismatic coefficient 
of 0 + 60. 

There is however a concomitant change in a parallel 
middle body found by solving for x at x = p, i.e.. 



Therefore the resulting change in prismatic coefficient is 
obtained by altering the length of parallel middle body and 
then proportionally expanding or contracting the entrance 
and run. Because of this procedure the method has the 
following disadvantages: 



simply reduces it by the factor 



1 ~ (0+60) 
1-0 





( 2 . 2 ) 
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1. There is no control over the length of the 
parallel middle body, i.e., <J> and p, cannot 
be varied independently. 

2. The procedure cannot be applied to reduce the 
fullness of a ship having no parallel middle 
body. 

3. Conversely, a ship cannot be increased in 
fullness without introducing parallel middle 
body. 

4. The prismatic coefficient of the entrance or 
run cannot be altered. - . 

5. The region where fullness is added cannot be 
controlled. That is, the maximum changes in 
fullness take place at the shoulders of the 
curve, i.e., point B. 

It is because of these numerous, severe restrictions that 
Lackenby sought to develop a more general technique of 
modification. 

2.2.2 Varying the Fullness of an Entrance or Run not 
Associated with Parallel Middle Body. 

In reference [ 2 ] , Lackenby concerned himself with 
providing a means by which to modify both C p , LCB and the 
length of parallel middle body in a controllable manner. 
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While many of these relationships are of importance, only 
those which apply more specifically to destroyer-type hull 
forms will be pursued in any detail. However, for the 
readers' convenience, the most generalized case of Lackenby's 
formulas is included as equations (2.10) through (2.13) in the 
last part of this section. For the in depth derivations the 
reader is referred to the original paper. Nevertheless, the 
following derivation is the foundation upon which all the 
subsequent relationships are based. 

In referring to figure 2.4 the various quantities have 
a meaning identical to those of the previous section. The 
only additional term requiring definition is k, defined 
mathematically as follows: 



fc'o’ 3 * 



The only other difference between figures 2.3 and 2.4 is that 
figure 2.4 represents a hull form not having parallel middle 
body, i.e., p = 0, and as a consequence the length of the 
entrance and run equals that of the half length of the ship. 

In referring to figure 2.4, it is recognized that in 
order to preserve the form of the parent at both the end of 
the ship, (x = 1) and the middle of the ship, (x = 0) an 
equation for 6x of the following form would suffice: 
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fix = cx(l - x) 



where c is an as yet to be determined constant. It may be 
seen in Appendix A that the relationship for fix is: 

fix = x (1 - x) (2.3) 

<p (l-2x) 




Clearly, the equation for fix is of the form of a second degree 
polynomial (parabola) whose values are 0 at x = 0 and x = 1, 
as desired. This relationship also shows that the amount by 
which any section in the parent is shifted is a function solely 
of the unchanged longitudinal position x and some as yet unknown 
value 54>. 
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At this point we must turn our attention and consider 
both the entrance and run concurrently if we are to lend some 
significance to the quantity If we desire to specify both 

C p and LCB, we are in essence placing a requirement on the 
area under the sectional area curve and its moment about some 
axis (say x = 0) . Since equation (2.3) applies independently 
to the entrance and run, it should be possible to select the 
6<J>'s of these respective regions such that when taken together, 
the ship has the values of C p and LCB desired. 

At this point we introduce the following quantities: 
z = the distance of the parent ship's LCB 
from midships, normalized by the half 
length, (positive forward) . 

6" z = the required shift in LCB to obtain that 
required for the daughter hull form 
(positive forward) . 

Prime (') - denotes derived forms. 

Subscripts - e = entrance or forward half-body, 
r = run or after half-body, 
t = a property describing the 
entire ship. 
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Referring to figure 2.5 above, we may interpret the require- 
ments on Cp and LCB mathematically as follows: 

C D = ♦ 

= <J> t + = (4> e + <f> r ) + ($<l> e + 6 <J> r ) (2.4) 

LCB =z _ _ 

z' - z + 6z 

Summing moments of areas 

z’ (<J> t + t ) = x e <J> 0 + k e 6<(> e “ fx r <l> r + h r 6<J> r ] 

z' = jp- ^ x e ^ e + ” f^^r + ( 2 « 5 ) 
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If equations (2.4) and (2.5) are rearranged, the expressions 
for <5<j> and 64 may be obtained. 



64 e = 


2 


’ {6<|, t (h r +2) + 


«2(* t + 


«* t )} 


(2.6) 


( W 


ii 

n 

*o 


2 


■ {«* t (h e -z) - 


6z(4 t + 


«* t )> 


(2.7) 


<h +V 



At this point the only variables which were not previously 
defined are h Q and h r . The exact expression for these 
variables is, with the appropriate subscript: 

h = 2^ + ii (a . 8) 

l-2x * (1 - 2xp 

While equation (2.8) contains 6<j>, the very thing for which it 
is being used to calculate, it has been stated [2 ] that the 
leading term along provides a very good approximation to h 
for "moderate" values of i.e., 

h = 2x T . 3k2 _ (2.9) 

1 - 2x 

Should it be desired to calculate 6<J> using equation (2.8), the 
solution will prove to be a quadratic which, while unwieldy, 
is certainly not unsolvable. The derivation of h may be found 
in Appendix A with the value of r defined as follows: 
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r 



3 




dy 



While the above equations with the derivations in the 
Appendix illustrate the underlying theory, the following 
expressions represent the most general form of Lackenby's 
work. 



S^e = a— +" b ~ {2 I 6 V B r + * )+6 * ( V 6 V 1+ V p e"V p r } 

( 2 . 10 ) 

6 *r = §" t ' B ^ 2 t<S<J> t (B e -z)-6i'((J, t -6(}, t )]-C e 6p e -C r 5p r } 
e r 

( 2 . 11 ) 

In the following expressions the items refer to the entrance 
or run as appropriate: 



«X - (X - X) {j|E + ^SZEL [64>-6p -|x=pfn 

The practical limits on 6<J> e and 6<J> r are: 

<5p (1-<J>) ± ^A[l - |£_] 

4 * = i=p 

A, B and C are calculated as follows: 

A = <J> (l-2x) - p(l-<jj) 

B = |{2x - 3k 2 - p(l-2x) } 

C = TTp {B(1 “ < < ,) " $ d-2x) } 



( 2 . 12 ) 



(2.13) 



(2.14) 

(2.15) 

(2.16) 
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It should be recognized that in the above equations the 
necessary section shift, 6x, is a function of the new values 
of Cp and LCB, and the properties of the original parent hull 
form. 

2.2.3 Modification of Lackenby's Method to Accommodate Hull 

Forms with Maximum Sections not at Midships . 

It should be realized that in all of the preceding 
developments it was assumed that the section of maximum area 
fell at midships. While this is certainly the case for a 
large class of vessels, it is virtually never true for contem- 
porary cruisers and destroyers. It is for this reason that a 
new set of equations was sought while still adhering to the 
basic philosophy of longitudinally shifting sections. 



HZ 
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Referring to figure 2.6, the definition of terms is, 
once again, consistant with the preceding section. There are, 
however, two very important changes. First, where z and fiz" 
were originally normalized by the ships half length, they are 
now normalized by twice that length or the length between 
perpendiculars, L. Second, x, the local longitudinal position 
of any section is normalized by the appropriate length of 
entrance or run L q or L r respectively. Also all values of 
x and z take as their origin the station of maximum sectional 
area.- 

The basic relationship for 6x is still of the form 
6x = cx(l-x) or 6x = ^ ' ^-^x) the same as equation 

(2.3) previously. However, equations (2.4) and (2.5) now 
become : 



<j> t = ^ + *H e ) + L r ^ r + (2.17) 

z' = z + 6z = {Lg(<fr e x e +5<|)e h e ) - I, r ( ^r 5F r +<S4, r h r )} 

L <J>. 

t (2.18) 
These two equations are solved simultaneously for 6<j> e and 
6<p in Appendix B, the results of which are listed below: 



e - .2, [j ~ ~ < l r*A- 1 I*eV ! ' L2 *;- 1 r h r ( Ve + Vr- i t;» 



= 



(2.19) 

5$ = 1 -L 2 <^^x -¥'L 2 *:-L h (L* +L r <f -L^) > 

L_h_+L_L_h_ eee rrr tee'ee rr t 

( 2 . 20 ) 
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I 

I It was these two equations along with equation (2.3) which 
proved to give very satisfactory results for several sample 
calculations. 

j 

2.2.4 A Method by which Constraints may be Placed on the 
Design Waterline 

In the preceding development the designer had no 
control over the shape of the design waterline. Because the 
longitudinal center of flotation LCF, was felt to be an 
important parameter in determining a ship's performance in 
a seaway, Moor [ 8 ] further developed the method of Lackenby 
to include control over the design waterline. It was with 
his revised method that Moor and his colleagues developed 
four distinctly different models with only their midships 
section identical. They were thus able to cut these four 
models in half to generate sixteen uniquely different hull 
forms. 

An interesting side light of this experiment was that 
the parent form used was that of a fast twin-screw currently 
in service whose maximum section was abaft midships. They, 
therefore, had to first swing the original area curve to 
place the maximum section at midships and then proceed with 
the new method of modification. Although having the maximum 
section at amidship may have proven to be more tractable for 
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the creation of the models, it was not necessary for the 
application of Lackenby's method. This fact was demonstrated 
in the previous section. 

The essence of Moor's method is that the sectional 

area curve and the design waterline are both altered in the 

sense of Lackenby to produce the desired values of C p , LCB, 

C w and LCF . At this point a new factor is introduced; the 

ratio between the sectional area ordinate and the design 

waterline ordinate is calculated for both the parent and 

derived hull forms. These ratios are then plotted as a 

function of ship length and it is from this curve that the 

longitudinal shift of sections is determined. Figures 2.7 

and 2.8 illustrate the sectional area and design waterline 

curves and the area/waterline ratio curve respectively. 

Referring to figure 2.8, to obtain the offsets for a 

particular section R, in the derived form, section R in the 
d p 

parent is used as a basis. The reason for selecting station 

R in the parent is that it is the closest section to R, with 
P a 

the same value of area/waterline ratio. The offsets of 

section R^ are then multiplied by the ratio of the beam 

coefficient in the derived form at section R, to that of the 

d 

parent at section R , i.e., B./B . These values may be seen 
P a P 

in figure 2.7. Additionally, if the maximum beam of the 
derived form is different from that of the parent, the 
offsets of section R^ also have to be multiplied by the ratio 
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of the maximum beam in the derived form to the maximum beam 



in the parent, i.e., ^max^pmax * Therefore, the equation 
for any offset b^ in the derived form is: 




^dmax 

b 

pmax 



( 2 . 21 ) 



It can be seen in figure 2.8 that there are regions of 
ambiguity. Such is the case where the derived curve lies 
below a minimum in the parent curve. It has been this 
author's experience confirming that in reference [ 8 ], that 
the regions which cannot be explicitly be defined may be 
faired in after defining the sections on either side. 

The one remaining undesirable characteristic occurs 
in regions where there is some form of keel rise, i.e., the 
fore foot or skeg region. In these areas, if the draft of 
the parent is proportionally altered to equal that of the 
derived form, there is a concomitant and undesirable change 
in the area of the section. It is to this matter which the 
next section is addressed. 
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2.2.5 A Method by which Constraints may be Placed on the 

Ship's Profile 

Since Moor's method proved capable of constraining 
both the sectional area curve and the design waterline, it 
was decided to extend the method to include the centerline 
profile of the ship. The actual mechanics require only the 
introduction of a local draft coefficient, (local draft/ 
maximum draft) into the denominator of the area/beam ratio. 
This new ratio, (Area/Beam/Draft) , is graphed and the sec- 
tional shifts determined from this graph. In determining 
the new offsets not only are the offsets of the parent 
modified transversely as described in the previous section, 
they are also altered in the vertical sense. This altera- 
tion is accomplished by using the water plane as a reference 
and moving the water lines below a distance proportional to 
the ratio of the derived form draft coefficient/parent draft 
coefficient. Also if there is a difference in the maximum 
draft of the derived form and parent, the waterlines are 
altered by this ratio as appropriate. 

As was mentioned in the background section of this 
chapter, section 2.1, for this modification technique to 
be mathematically rigorous the sections of maximum area, 
maximum design waterline breadth and maximum draft must be 
coincident. If this is not the case, the actual areas of 
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the sections generated will be consistently different by a 

very small amount from what is desired. From the few 

examples this author has worked, it is estimated the 

difference in the value of C obtained and that desired is 

P 

on the order of 1%. The explanation of this is seen in 
Appendix C. 
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3. Mathematical Representation of the Lines of a Ship 



3.1 Background 

It was established in section 2.1 that before the 
naval architect attempts to actually draw the lines of a new 
ship, he must have a "firm" description of the new design as 
represented by the various coefficients and curves of form. 
Examples of these, as cited previously, are: the sectional 
area curve, and hence Cp and LCB, the design waterline curve, 
(C w and LCF) , the principle dimensions and perhaps specific 
information about the geometry of the midship section or 
stern region. These characteristics should represent what the 
designer feels is the "optimum" solution to his set of re- 
quirements. The naval architect now has to create one or 
more, of a possible infinity of, design candidates which ful- 
fill his descriptive coefficients. 

The traditional method for drawing the various lines 
of the ship is with the use of long, continuous strips of 
wood, metal, or more recently plastic, held in the desired 
position by weights. These tools are called splines and 
ducks respectively. The curves produced by this method 
were continuous but often times contained unwanted waviness. 
Removing these unwanted undulations, while preserving the 



40 



desired character of the line, is a process referred to as 
fairing. This topic, and the implications of placing a 
mathematical interpretation on it, are discussed in the 
next chapter. Not only did the naval architect have to 
generate smooth curves which pleased him visually, there also 
had to be a consistancy in location of the surface points 
when observed from the different views. This is sometimes 
referred to as cross-fairing and is also addressed in the 
final chapter. It is the fairing, and cross-fairing, which 
represents a very large part of the manual design effort. 

It was recognized long ago, that if the ship design 
process was to be automated to any degree, a technique to 
represent the lines of the ship mathematically would have 
to be developed. This is especially true today where much of 
the work is to be done by high speed digital computers. Not 
only must the designer/programmer provide the mathematical 
algorithms for representing the ship's lines, he must also 
provide the logic necessary for the computer to duplicate the 
heretofore trial and error methods of the draftsman. The 
alternative to programming the logic however, would be to 
give the system an interactive man-machine interface at the 
decision points. It is, however, the mathematical repre- 
sentation of these ships' lines to which this chapter is 
devoted. 
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3.2 Development 



3.2.1 Derivation of the Spline Cubic Equation using a 

Variational Approach 

There are essentially two different methods by which 
one may arrive at a mathematical representation for ships 
lines. 

1. Select some mathematical function with several 
unspecified parameters whose values may be deter- 
mined by some accuracy criteria and boundary 
conditions. Typical of this approach is the use 
of a polynomial and a least squares fit criterion. 

2. Choose some smoothness and closeness of fit 
criteria such that, when taken together with the 
boundary condition, the function and parameters are 
determined. 

It is this second method, based on a variational smoothness 
criterion, that will be developed in this chapter. 

In general these variational methods involve the 
minimization of the integral of some linear combination of 
the squares of the various derivatives of the function 
sought. In the case where the equation of the flexible 
spline is sought, the "smoothness" criterion is taken to be 
the minimization of the strain energy in the spline. 
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